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We provide upper bound on the maximal rate at which irreversible quantum dynamics can gen¬ 
erate entanglement in a bipartite system. The generator of irreversible dynamics consists of a 
Hamiltonian and dissipative terms in Lindblad form. The relative entropy of entanglement is cho¬ 
sen as a measure of entanglement in an ancilla-free system. We provide an upper bound on the 
entangling rate which has a lnarithmic dependence on a dimension of a smaller system in a bipartite 
cut. We also investigate the rate of change of quantum mutual information in an ancilla-assisted 
system and provide an upper bound independent of dimension of ancillas. 

PACS numbers: 03.67.Bg, 03.65.Aa, 03.67.Hk, 03.67Mn 


INTRODUCTION 

The problem addressed in this paper is, given some 
Hamiltonian and dissipative interactions between two or 
more subsystems, what is the maximal rate at which an 
entanglement can be generated in time. The problem of 
upper bounding the entangling rate is called small in¬ 
cremental entangling. See the next section for precise 
definition. 

The same question for closed bipartite system evolv¬ 
ing under a unitary dynamics was answered by Bravyi in 
[1] and by Acoleyen et al. in [2] and Audenaert [3] for 
the general case in the presence of ancillas. Let us say 
that two parties, Alice and Bob, have access to systems 
A and B respectively together with ancilla systems a and 
b respectively. The system starts in a pure state \'&) aAB b 
and evolves according to a Hamiltonian Hab that acts 
only on systems A and B. Since the state of a system 

I ^WtaABb = exp {itH AB } \^) aABb sta y s P ure > one ma Y 
calculate the entanglement entropy as a measure of en¬ 
tanglement in a bipartite cut between Alice’s aA and 
Bob’s systems Bb, E{t) = — Tr(p aj 4 (t) In p aA (t)), where 
PaA{t) = Tret, |’F(t)} {'&(t)\ aABb . The entangling rate is 
a time derivative of the entanglement entropy at time 
t = 0, T(T, H) = dE{p{t))/dt\t=o- The small incremental 
entangling problem aims at finding an upper bound on 
the maximal entangling rate T(iJ) = sup^ T(T, H) 
that is independent of dimensions of ancillas and the ini¬ 
tial state |^ r ) aj 4Bb i The problem of maximizing the en¬ 
tangling rate of a bipartite system in the presence of an¬ 
cillas evolving under unitary dynamics was studied by 
many authors. 

In a case when A and B are qubits, Childs et al [4] 
gave upper bounds for an entangling rate and showed 
that they are independent of ancillas a and b. Wang and 
Sanders [5] considered systems A, B and ancillas of arbi¬ 
trary dimensions and proved that the entangling rate is 
upper bounded by T(H) < j3 ~ 1.9123, for a self-inverse 
product Hamiltonian. Bravyi [1] proved that in a gen¬ 
eral case with no ancillas the entangling rate is bounded 


by T{H) < c(d)\\H\\ lnd, where d = rain(d A ,dB) is the 
smallest dimension of the interacting subsystems and c is 
a constant close to 1. For an arbitrary bipartite Hamil¬ 
tonian Bennett et al [6] proved that the upper bound on 
the ancilla-assisted entanglement is independent of the 
ancilla dimensions T(H) < cd 4 ||iij|, where c does not de¬ 
pend on a or b. The bound was improved by Lieb and 
Vershynina [7] providing an upper bound T(H) < 4||7J||d 
for an arbitrary Hamiltonian in ancilla-assisted system. 
Finally the question was answered by Acoleyen et al. in 
[2] arriving at T(i7) < 18||.ffj| lnd. Few months later an 
independent proof was presented by Audenaert [3] that 
gives an upper bound T(H) < 8||.ffj| lnd. The bound with 
lnarithmic dependence is optimal, since one can find a 
particular Hamiltonian H A b for which they there is an 
equality. See Marien et al. [8] for a detailed review on 
the entangling rates for bipartite closed systems. 

In this paper we go further by considering open systems 
evolving under irreversible dynamics and aim to provide 
an upper bound on the entangling rate of a bipartite sys¬ 
tem. The behavior of entanglement under non-unitary 
evolution was first studied by Rajagopal and Rendell [9] 
and Zyczkowski et al [10]. Rajagopal and Rendell studied 
a dynamic evolution of a pair of initially entangled har¬ 
monic oscillators in the presence of local environments. 
Zyczkowski et al considered a bipartite system consist¬ 
ing of two spin-1/2 particles, one of which is subjected 
to periodic actions of a quantum channel representing 
the interaction with environment, and the whole system 
undergoes a sequence of global unitary evolutions. Both 
papers found that the entanglement may exhibit revivals 
in time and that it might vanish at finite times. For fur¬ 
ther research on revival of entanglement and its collapse 
see, for example, [11]-[13]. Another motivation to study 
entanglement rate problem in open system arises from 
investigation of production and robustness of entangle¬ 
ment in open systems, see, for example, [14], [15]. See 
[16] for a comprehensive review of the study of dynamics 
of entanglement in open systems. 

To investigate the entangling rate in an open system 
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we can not consider an entropy of entanglement as an 
entanglement measure, since for a system starting in a 
pure state and evolving under non-unitary evolution a 
time-dependent state may no longer stay pure. 

To quantify the entanglement in ancilla-free system we 
take the relative entropy of entanglement, which reduces 
to the entropy of entanglement on pure states. We derive 
an upper bound on the entangling rate that has a lnarith- 
mic dependence on the dimension of a smallest system, 
as it was in the case of a unitary dynamics. These result 
rises a question of finding a class of entanglement mea¬ 
sures for which the same holds true together with another 
question whether the same bound can be derived in the 
ancilla-assisted case. 

In ancilla-assisted system we consider quantum mutual 
informal and provide an upper bound on its derivative 
that is independent of ancillas. 

The paper is organized as follows: in Section 2 we dis¬ 
cuss a state evolution in an open system, talk about en¬ 
tanglement measures, describe the small incremental en¬ 
tangling problem in a closed system and conjecture the 
problem for the open system. In Section 3 we provide an 
upper bound on the entangling rate for a relative entropy 
of entanglement in ancilla-free system that has lnarith- 
mic dependence on the smaller system. In Section 4 we 
derive an upper bound on the derivative of quantum mu¬ 
tual information in ancilla-assisted system. 


PRELIMINARIES 

Suppose that two parties, say Alice and Bob, have con¬ 
trol over systems A and B. In ancilla-assisted entangling 
both parties have access to additional subsystems, called 
local ancillas, i.e. Alice is in control of two systems A 
and a and Bob is in control of B and b. 

Alice and Bob start with a pure state p a ABb( 0) = 
Ilk} (^\ aABb - A time dependent joint state of Alice and 
Bob p(t) satisfies the following differential equation 

j t p(t) = C AB {p(t)). ( 1 ) 

The generator of the dynamics Lab acts only on systems 
A and B, and consists of both Hamiltonian part H A b and 
a dissipative part in Lindblad form with operators Lab 
and takes the following form [17] 

L A B(p)=-ilH A B,p\ + J2( L * B P( L * B y ( 2 ) 



We are going to consider entanglement measure E(-) 
that satisfy the following assumptions: 

1. E vanishes on product states 

2. E is invariant under local unitary operations 

3. E can not increase under LOCC operations 

4. for any p we have E{p ® |$d) (4>d|) = E(p) + 
£(l$d)), where |$ d ) = ^ Ej=i \j) 1 \j) 2 is the 
maximally entangled state shared between two par¬ 
ties. 

There are several entanglement measures that sat¬ 
isfy these conditions. In this paper we will focus on 
entanglement measures that reduce to the entropy of 
entanglement on pure states 

5. for p = \(j>) {(j> | 12 , we have E(p) = —Tri(pi lnpi) = 
-Tr 2 (p 2 hr pi)- 


In case of a unitary dynamics, i.e. L a = 0 for all a , 
the total change of entanglement though time is no more 
than 21nii with d = min{dA,ds}; see [6], [8]. Property 
4 of the entanglement measure plays a significant role in 
the bound of the total change of entanglement. In the 
case of irreversible dynamics the total change of entan¬ 
glement depends on the measure of entanglement and it 
is unknown for many measures. For example, a relative 
entropy of entanglement (see definition in (4)) gives an 
trivial bound of 41nZI with D = min{d a A, d B b}, and neg¬ 
ativity A f gives a linear bound in d a A or d B b , see Tiersch 
et al. [18]. 

A unitary case gave motivation to derive a bound on an 
instantaneous rate of change of the entanglement. This 
problem is formulated using the entangling rate. The 
case of a unitary dynamics entangling rate is defined as 


r (*,h) = 


dE{p{t)) 


dt 


i=0 


For a system evolving according to a unitary evolution 
the following theorem was conjectured in [1] and was first 
proved in [2]. 

Theorem 1. [2] Small incremental entangling for 
unitary dynamics. Denote d = min{ dA,d B }. For a 
system evolving according to a Hamiltonian H AB there 
exists a constant c such that the entangling rate is 
bounded above by 


Here H AB , {L^ B }a are bounded linear operators and 
H AB is also Hermitian. We will be dropping the (sub) su¬ 
perscription notation AB throughout the paper for easier 
notation when it is clear from the context which systems 
an operator acts on. 


r(<F,Ff) < c||J?||lnd, 

where c is independent of the dimensions of systems A 
and B, ancillas a, b, the Hamiltonian H AB and the initial 
state IT). 
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The upper bound on the maximal entangling rate was 
investigated in many papers [1], [4]-[7] and was first 
proved in [2] with a different proof provided in [3], see 
also [8] for detailed review of this case. 

We ask the same question for a system evolving under 
an irreversible dynamics (1). In general case an entan¬ 
glement measure may not be differentiable. Define the 
entangling rate for time interval At > 0 as 


WAA t) 


E(p(t)) - E(p( 0)) 
At 


In case when the entanglement measure is differentiable, 
the entangling rate is defined as 


r(*,£) 


dE(p(t)) 


The relative entropy of entanglement [19] of a state 
p(t) is given by 

D(p{t)) : = mmD(p(t)\\a) (4) 

crsep 

= minTrfp(t) In p(t) — p{t) In ex'), 

crsep \ / 

where a = J2j a j a A(j) ® cr B (j) with J2j a j = 1 is a 
separable state between Alice and Bob’s systems. For 
pure states p := p(0) = |’F) ('Ll the relative entropy of 
entanglement is an entropy of entanglement of this state 

[ 19 ]- 

In this section we prove the following theorem. 

Theorem 2. For any e > 0 there exists 8 > 0 such that 
for any At < 8 the entangling rate for the relative entropy 
of entanglement has the following upper bound 


Similar to a system evolved under unitary evolution we 
are looking for the upper bound on entangling rate that 
is independent of ancillas and has a proper dependence 
on dimension of systems 

Question 1. Small incremental entangling. De¬ 
note d = min {d A ,d B }- For which entanglement mea¬ 
sures there exists a constant c and a non-negative non¬ 
decreasing function /(•) such that for any e > 0 there 
exists 8 > 0 such that for any At < 8 the entangling rate 
is bounded above by 


T*(*, A At) < 4(||tf|| + 86 Y, ll^ll 2 ) Ind + e, 

Ot 

where d = min(d J 4 , d B )- 

First, let us look closely on the relative entropy of en¬ 
tanglement for pure states. Let state |\H) have Schmidt 
decomposition 

d 

\^) = ^VP^\ ( t ) n) A \^n) B - ( 5 ) 

n =1 


r(AAAt) <c\\C\\f(d) + e, (3) 

where c is independent of the dimensions of systems A, 
B, ancillas a, b, the generator C and the initial state 

\^)aABb' 

Note that the intuitive argument to consider the open 
system dynamics on system AB as a part of the unitary 
dynamics on the bigger system ABE doesn’t work in 
this case, because there is no bound of the type ||£ab|| > 
C\\F ABE \\ with C > 0 being independent of the generator 
C A b and the induced Hamiltonian F A be- The bound of 
this type does not exist because in a case when C AB = 0 
(.H = I, L a = 0) and F AB e = Iabe constant C must be 
zero. 

We answer Small Incremental Entangling question for 
the relative entropy of entanglement in an ancilla-free 
case and for the quantum mutual information in ancilla- 
assisted system. 

RELATIVE ENTROPY OF ENTANGLEMENT IN 
ANCILLA FREE SYSTEM 

Assume throughout this section that there are no an¬ 
cillas, d a = db = 1, and that d := d B < d A . For states 
and operators that act on both systems A and B we will 
drop the subscript notation AB. 


Then the relative entropy of entanglement is achieved by 
a state [19] 

d 

(JO = ^Pn \4>n) {f>n\ ® \4>n) (Ar| • (6) 

71—1 

Proposition 1. For states p = |4/) (T| and do defined 
in (5) and (6) respectively, there exists a mixed state p 
such that 

= G+ (1 ~ 

Proof. The statement of the proposition is equivalent to 
inequality 

dcr 0 > p = |'F) (’Ll, (7) 

which in turn is equivalent to inequality 

Z ■= d^Pn | (fni’n) y/PFPk \(f>n4>n) {4>ki>k\ > 0. 

n n,k 

The last inequality is equivalent to a statement that 
for any state, which can be written as |fl) = 
Enfc V ank IV’fc)) the expectation value of the observ¬ 
able Z is positive, i.e. 

(H| Z |fl) = dy^ y p n a nn — (£ y/Pn^nn^ 0* (8) 
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Since a function f{x) = x 2 is convex, we have that for 
any x n ,n = 1, 

n—1 n 


following form, for any 0 < p we have 

Tr(L t [Lp,lnp]) = -Tr(£ f [.L(pp),hip]) 
P 

=: ?-Tr(tf[LX,lnY]), 


Taking x n = y/p n oinn in the above inequality we arrive 
at inequality (8), which in turn gives (7), proving the 
proposition. 

□ 


Now we are ready to prove Theorem 2. 


Proof. The relative entropy of entanglement of a state 
p(t = 0) at time zero is the relative entropy of that state 
D{p) = D(p\\ao) = E(\ h), where <r 0 is defined in (6). For 
any time t we have that D(p(t )) < D(p(t)\\cro). Therefore 
for any e > 0 there exists <5 > 0 such that for any At < S 


r R (%£,At)<±D(p{t)\\v 0 ) 


+ e. 
t=o 


(9) 


The derivative of the relative entropy D(p(t) | |cr 0 ) can 
be calculated as follows 


4 ; D (p(t)\Wo) 
dt t =o 

= Tr(p(t)| t= olnp — p(t)| t= o In cr 0 ) 

= —*Tr([iJ, p] In p) + iTr([H, p] In cj 0 ) 

+ 5Z Tr ( L “^ L a ln ^~ \{ L a L «,p} Inp) 

Ot 

- 5Z Tr ( ictpL “ lncro ) +5I Tr (^ L “ ia ’^ lncro ) 

a. a. 


iTr^H\pp,\na 0 fj 

(10) 

- Tr (pp) ’ in 

' ot 

(11) 

- Tl '( L *[(,PP) L li lllCTo]) 



ot 


- 5Z Tr ( L “[ ia ^, ln ^)' 


Here for the first equality we used that Tr(/j(f)|t = o) = 
Tr(£(p)) = 0, which can be seen from the expression of 
the Lindbladian generator (2) and the cyclicity of a trace. 
For the third equality we denoted p := 1/d. Note that 
the third term in the last equality is a complex conjugate 
of the second one 


Tr(pt[PX, In Y]) = Tr(P[XPt, InF]), 

where X := pp, Y := <7 q = pp + (1 — p)p according to 
Proposition f. The fourth term can be brought in the 


where 0 < X < Y, TrX = p, TrF = 1. Therefore 


—D{p(t)\\a 0 ) 


<-\Tr(H[X,lnY])\ 
t=o P 


+ 2^|Tb(4[T a X, 1 ny])|, 


where TrX = p , TtY = 1 and 0 < X < Y. 

Consider two terms (10) and (11) separately. 

1. Term (10) is linear in H , so we may pull out the 
norm ||F7||, then the term becomes exactly the one that 
arises in a unitary case 


Tr(P[X,lnF]) (12) 

for ||P|| = 1, TrX = p, Tr Y = 1 and 0 < X < Y. 
Acording to [3] the following upper bound holds 

Tr(P>p,ln<7 0 ]) < -2[p\np+ (1 — p)ln(l -p))||Pj|. 

For any p < 1/2 one has — ^plnp + (1 — p) ln(l — P)) < 
—2plnp and therefore 

Tr(P[p,In(7 0 ]) < -41np||P|| = 41nd||P||. (13) 

2. Pulling out the norms of L a term (11) becomes 

Tr(pt [LX, In F]) = \\L\\ 2 Tr(D[LX,lnY]), 

with ||L|| = 1, TrX = p , TrF = 1 and 0 < X < Y. 

Any operator can be written as a sum of hermitian and 
anti-hermitian parts, let U = 2P — I + i(2R — I) with 
0 < P, R < I. Denote Q := LX , then 

|Tr(pt[Q,l n F])| < \2Tr{P[Q,lnY]) + 2iTr(R[Q,lnY])\ 
<4|Tr(P[Q,lnF])|, 


where 0 < P < I denotes either P or R. 

Expression |Tr(P[Q, lnF])| is similar to the expression 
of the entangling rate for a unitary dynamics (12), which 
allows us to follow the argument presented in [2] to ob¬ 
tain an upper bound on this type of expressions. For 
completeness sake we present a full analysis here, which 
globally resembles an argument in [2], but since Q = LX 
is not even Hermitian some details need to be different 
in order to derive an upper bound. In the future we will 
drop tilde above L in expression for Q. 

Let yj denote the eigenvalues of Y in decreasing or¬ 
der 1 > yi > y 2 > ■■■ > Vn > H with corresponding 
eigenstate \4>j). Denote the matrix elements of Q in this 
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eigenbasis as Qij = Q \4>j), similarly define the ma¬ 
trix elements of P as Py = (<j>i \ P\(j>j). Then 


|Tr(P[Q,lnF])| = I £ In ^(QyP* 
1 


QjiPij) ■ (14) 


Group the eigenvalues of Y in the following intervals 


1 > 2/ii > P 1 < ii < ni 

P>Vi 2 > P 2 m < Vi2 < n-2 

(15) 

P k ~ l > Vi k >P k n k <i k < N. 

Some of the intervals can be empty, in this case then 
n k ~ i = rife. Now rearrange the sum (14) in the following 
way, for A j € (n k -i,n k ) 

£ = £ + £ +■■• + £ 


l<il<jl<n 2 

ni<i 2 <j 2 <n 3 

nk-l<ik<jk<N 

(16) 

£ 

- £ - 

£ 

ni<i 2 <j 2 <n 2 

n 2 <i 3 <j 3 <n 3 

"li -1 <ik<jk<n k 

(17) 

( E + 

£ +■■■ + £ 

). (18) 


Al,Afc>2 A2,Afc>3 Afe_2,Afc 


The first two lines contain sums that run over the one 
or two consequent intervals in (15). The sum on the last 
line runs over pairs that are separated by at least one 
interval in (15). 

The sums in the first line (16) can be calculated in the 
following way. Let us bound, for example, the first sum 

712 712 

I £ £ \n^-(Q ij P ji -Q ji P ij )\ = \Tr(P[Q,lnY])\ 

7=1 j=i+ 1 ^ 

< 2|Tr(P[X, In ^-])| + 4|Tr(P[TX, In (19) 

2/771771 2/777771 

Here we wrote out operator L = 2P — I + i(2R — I) as a 
sum of hermitian and anti-hermitian parts and denoted a 
Hermitian operator T as either operator P or R. We also 
added a zero term by scaling Y. Tilde above an operator 
denote a restriction of this operator on a space spanned 
by eigenstates \(pj) for j = 1, ...,ri 2 . 

The first term in (19) is bounded by 

2|Tr(P[X,ln^£])|< ||[X,ln^£]|| 1 
2/777777 Umin 

< ln£^||X||i 

2/777777 

< 2(pi + p 2 )ln(l/p). 

Here for the first inequality we used the fact that for any 
Hermitian operator A 

M||i = max |Tr(PA)|=2 max |Tr(PH)|. (20) 

||P||<1 0<P<I 


For the second inequality we used Kittaneh inequality 
[20] for a positive compact operator A and a unitarily 
invariant norm |||[H,X]||| < ||H|| |||Jf|||. On the last in¬ 
equality we used p 2 < yi/yj < 1/p 2 and ||^T||! = TrX = 
Pi + P 2 , where TrX = J2 k Pk = P- 

The second term in (19) can be bounded similarly 


4|Tr(P[fX,ln 




_ Y _ Y 

< 4|Tr [PTX In--)| + 4|Tr(TXPln --)| 

Vmin Vmin 

4|| In -7—— 1| (\TY(W PTX)\ + \Tr(PWTX)\) 

Um.i.n. ' ' 



Vmin 


<8\n Vrnax \\ 


Vmin 

Here 

for the 

In tX 

Umi 

-/n in s£n 


ity (20) for the second inequality, since ||HLPT|| < 1 and 
\\PWT || < 1. 

Therefore line (16) can be bounded above by 


fc-i 


18 (pi +p k + 2 £pj) ln(l/p) < 36pln(l/p). (21) 


i=2 


Line (17) can be bounded in a similar fashion. The 
first term, for example, gets a bound 


Pr, 


n 2 n 2 

| £ £ In^iQijPji - QjiPij)\ < 5Tr(X) In ■ 

, , 2// Vmin 

7=77i J — 7+1 


< 5p 2 ln(l/p). 

Therefore the whole second line (17) can be bounded by 

5(p 2 + ... +pfc)ln(l/p) < 5pln(l/p). (22) 

To bound the last line (18) define Z = Y~ 1 / 2 QY~ 1 ^ 2 . 
Then the last line (18) reads 

|£ . In yiyl /2 y 1/2 (Z ijPji - ZjiPij )| 

z — l< 3 Vj 


< 


(E,<>| ypyy^z-, 

*(t<> 

+ (£.<,■ 


1/2 

1/2 

1/2 

1/2 


< 


+ 


N 1 , 2 N 

P 1/2 ln(l/p) ( £ 2 liZijiZ^jYJ ( £ y, 


\Pij\ 


i,j =1 


*>J=1 


N 1/2 N 

p 1/2 ln(l/p)( £ y,Z,,[ZA),^ ( £ P 


i\Pji 


1/2 

(23) 

1/2 


*J =1 


7=1 
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Here for the first inequality we used Cauchy-Schwarz 
inequality. For the second inequality we used that 
x 1 / 2 ln(l/x) < p 1 / 2 ln(l/p) if x < p < 1/e 2 , with 
x = Vj/Vi- 

Since 0 < P < /, we have that 0 < P 2 < / and 
therefore for every i we have 1 > ( P 2 )u = Y/j PijPji = 
J2j I P-ij 1 2 ■ Using /T, '(/, = Try = 1 we continue calcula¬ 
tions 

(23) < p 1/2 1 n( 1 /p)( X r(yZZ t ) 1 / 2 +Tr(YZ^Z) 1 / 2 \ 

(24) 

Considering each trace separately, recall that Z = 
Y-^LXY- 1 / 2 . Then 

Tr (YZZ f ) = Tr (PLXY^X). 

Denote R := XY^X. Since 0 < L^L < I, 0 < 
X 1 / 2 Y~ 1 X 1 / 2 < /, we have that 0 < R < X and 
RP' 2 (I - lM)R}' 2 > 0. Therefore 

TrfYZZ^) = Tr (I3LR) < Tr(P) < Tr(A) = p. 

The second term in (24) can be estimated in the following 
way 

Tr (YZ^Z) = Tr (XL^Y~ 1 LX) =: Tr(AU), 
where A := Y~ 1 / 2 LX. Therefore 
Tr (YZ f Z) = £ |A,(A)| 2 = £ |A .(X 1 ' 2 AX- 1 ' 2 )] 2 

3 3 

= Y |A j (X 1/2 Y~ 1/2 LX 1/2 )\ 2 =: Tr(ptp), 
j 

where A j(A) is j-th eigenvalue of operator A and B := 
X 1 / 2 Y~ 1 / 2 LX 1 / 2 . Here we used the fact that A and 
P 1 / 2 AP~ U 2 have same eigenvalues for any P since they 
have the same characteristic polynomial. Therefore writ¬ 
ing out pip we obtain 

Tr -(YZ f Z) = Tr(y- 1/2 Xy- 1 / 2 PA:P t ) < Tr (LXL f ) 

< Tr(X) = p, 

here we used that LXLf = (X 1 / 2 pi)i(X 1 / 2 Pi) > 0. 
Thus line (18) is bounded above by 

(18) < 2pln(l/p). (25) 

Combining bounds (21)-(25) all lines (16)-(18) in sum 
(14) are bounded above by the following expression 

\Tr(tf[LX,lnY))\ < 172||P|| 2 pln(l/p), (26) 

where p = 1/d. Using this bound together with (13), 
entangling rate (9) is therefore bounded above by 

r*(tf,£,Ai) <4(||P||+86^||P a || 2 ) In d + e, 

ot 

where d = min(d.A, d B )- □ 


Note that in the ancilla-free system if all Lindblad 
operators L a vanish, i.e. in the case of unitary evolu¬ 
tion, Theorem 2 reduces to Theorem 1 with constant 
c = 4. This constant comes from the use of a gen¬ 
eral bound provided by Acoleyen [2], which is called a 
small incremental mixing: for two states p\ and p 2 , any 
0 < p < 1, and any Hamiltonian H the following bound 
holds Tr(P[ppi, ln(ppi + (1 - p)p 2 )}) < —2{p\rrp + (1 — 

p) ln(l — p)^j. We used this bound in (13) with p = 1/d to 

derive Tr (H[pab, ln<T 0 ]) < 41nd||P||. If we consider rel¬ 
ative entropy of entanglement in the presence of ancillas 
we would have 1/p = min{d a A, d B b}, which depends on 
the dimension of ancillas. In ancilla assisted case small 
incremental mixing problem proves small incremental en¬ 
tangling problem with four times bigger constant, i.e. 
c = 8. Our ’’improvement” of the constant came from 
the fact that we consider ancilla-free system. The best 
known constant for the ancilla-free system was given by 
Bravyi [1], which is close to one. We could have used 
this bound instead of (13), but we chose not to since no 
precise expression of this constant was given. 

QUANTUM MUTUAL INFORMATION 

As a measure of correlations between Alice’s and Bob’s 
systems one may calculate quantum mutual information 
between these two systems. The quantum mutual infor¬ 
mation of a state p a ABb in a bipartite cut Alice—Bob is 
as follows: 

I{aA\ Bb) p = S(paA) + S(pBb) ~ S(p a ABb) 

= D(paABb\\p a A ® PBb )• 

Since mutual information is differentiable one may inves¬ 
tigate the rate of change of quantum mutual information 
for system evolving according to the irreversible evolu¬ 
tion (1). 

Theorem 3. For a system starting in pure state p a ABb = 
|4') (^\aABb an d evolving according to evolution (1) the 
following inequality holds 

—I(aA] Bb) p (t) 

ai t=o 

< 4(2||P|| + 129^ ||P Q || 2 )(lndA + \nd B ). 

Ot 

Proof. Taking a time derivative of a quantum mutual in¬ 
formation we obtain 

fjfI(aA', Bb) p{t) 

- Tr(p aA Bb(t)\t=o In PBb) 

+ Tr (p a ABb ( t) 1 1= 0 In PaABb)- 


= - Tr (paABb(t)\t=o In Pa a) 

t.= 
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Writing out the dynamics of a state according to (1), we 
arrive at 


Bb) p{t) 

ai t=o 

Ir 

= Tr(iH AB [p a ABM{paA ® 3 — )]) 

as 

~ \ ^2 Tr(Ll[L a p aAB ,ln(p aA ® ^-)]) 

a 

+ ^ Tr ( L « [PaAS-^, lll(p aA <gj ^)]) 
a 

+ Tr(iiJ y 4B[pABb,ln(-4- ® psb)]) 

d A 

~\Yh T r( L il L a PABb, In {j- ® PBb)]) 

O' 

+ 0 ^)]) 

a 

- y~lTr(L^[L a p aj 4B b , h\p aABb \). 

a 


According to Lemma 1 in [1] there exist states p a AB and 
<JABb such that 


„ Ib 

PaA ® -T- 

d B 


j2 PaAB T (1 ,2 )paAB 

d B d B 


and 


-r- ® PBb = -yrpABb + (1 “ 35 ~) a ABb- 

d A d\ d z A 

Following arguments discussed in previous section lead¬ 
ing to (26) we find that 


-iTr{H AB [p aAB ,ln(p aA ® 3 —)]) < 8 In d B ||-f7||, 

dB 


\Tx{Ll\L a p aAB Mp aA })\ < 172||L Q || 2 ln(4) 


and 


-iTr(H AB [p AB b, 111(3- ® PBb)]) < 8 lnd^H-HH, 
d A 


\Tr(Ll[L a p ABb ,\np B b})\ < 172||L a || 2 ln(4). 

Therefore the derivative of the mutual information is up¬ 
per bounded by a quantity independent of dimension of 
ancillas 

4 1(aA;Bb) p{t] 

ai t—0 

<4(2||i7||+129^||L Q || 2 )(lnd A + lnd B ). 


□ 


CONCLUSION 


In this paper we investigated a question of the exis¬ 
tence of an upper bound on the entangling rate of an 
open bipartite system that evolves under an irreversible 
dynamics. We considered two cases. First we considered 
an ancilla-free system and took a relative entropy of en¬ 
tanglement as an entanglement measure. We proved an 
upper bound on the entangling rate similar to the one 
in a unitary dynamic case. Having this bound we at the 
same time have obtained an upper bound on entangling 
rates for all entanglement measures that are majored by 
relative entropy of entanglement and reduce to the en¬ 
tanglement entropy on pure states, like entanglement of 
distillation [21]. Second, we discussed the rate of change 
of quantum mutual information in ancilla-assisted system 
and provided an upper bound independent of ancillas. 

A question of finding an upper bound on the entan¬ 
glement rate for any entanglement measure, or a class of 
entanglement measured, remain open. The small entan¬ 
gling rate problem was applied to the area law problem in 
the closed system [2], [8]. It remains to be seen whether 
one may find a similar behavior in an open system. 
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